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ABSTRACT 

This investigation focuses on the development of multibody analytical models to predict the dy- 
namic response, aeroelastic stability, and blade loading of a soft-inplane tiltrotor wind-tunnel 
model. Comprehensive rotorcraft-based multibody analyses enable modeling of the rotor system 
to a high level of detail such that complex mechanics and nonlinear effects associated with control 
system geometry and joint deadband may be considered. The influence of these and other non- 
linear effects on the aeromechanical behavior of the tiltrotor model are examined. A parametric 
study of the design parameters which may have influence on the aeromechanics of the soft-inplane 
rotor system are also included in this investigation. 


APPROACH 

The objective of this investigation is to develop 
and refine multibody analytical models to predict 
the dynamic response, aeroelastic stability, and blade 
loading of a soft-inplane tiltrotor wind-tunnel model. 
Comprehensive rotorcraft-based multibody analyses 
enable modeling of the rotor system to a high level 
of detail such that complex mechanics and nonlinear 
effects associated with control system geometry and 
joint deadband may be considered. The influence of 
these and other nonlinear effects on the aeromechan- 
ical behavior of the tiltrotor model will be examined. 
A study of the design parameters which have influ- 
ence on the aeromechanics of the soft-inplane rotor 
system has been addressed in a previous work, and 
its prosecution is part of this investigation. This work 
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also focuses on forward-flight configuration analysis, 
investigating stability results obtained from past ex- 
primental campaigns. This research is being per- 
formed as a cooperative agreement between the U.S. 
Army, NASA Langley Research Center, and the Uni- 
versity Politecnico di Milano. 

A new four-bladecl semi-articulated soft-inplane 
(SASIP) rotor system, designed as a candidate for 
future heavy-lift rotorcraft, was tested at model scale 
on the Wing and Rotor Aeroelastic Testing System 
(WRATS), a l:5-scale aeroelastic wind-tunnel model 
based on the V-22. Previous investigations involved 
a three-blade soft-inplane hub and mainly addressed 
the stability properties of this configuration (1, 2). 
Soft-inplane design implies reduced hub loads and 
should allow significant weight reduction in large 
scale tiltrotor aircraft. However, there is significant 
potential for reduced whirl-flutter stability margins 
in comparison to the currently exploited stiff-inplane 
configurations, so extensive aeromechanical investiga- 
tion is required to determine the feasibility and the 




Figure 1: MBDyn multibody model of the three- 
bladed stiff-inplane tiltrotor system. 

requirements of soft-inplane design. 

The experimental part of this investigation in- 
cluded a hover test with the model in helicopter mode 
subject to ground resonance conditions, and a for- 
ward flight test with the model in airplane mode 
subject to whirl-flutter conditions. A three-bladed 
stiff-inplane gimballed rotor system, used in several 
previous experiments, was examined under the same 
conditions as the four-bladed soft-inplane hub to pro- 
vide a baseline for comparison. 

Detailed analytical models of the SASIP tiltro- 
tor have been developed using two multibody rotor 
codes, one known as MBDyn (Ref. 3) and one known 
as DYMORE (Ref. 4) . The two codes have similar ca- 
pabilities, but it is desirable to compare their results 
(with models created by two different researchers at 
two different institutions) as a test of robustness for 
the multibody approach. The multibody analyses in- 
clude dynamic models for parts of the rotor system 
which are often not considered in classical rotor anal- 
yses, such as the hydraulic actuator control system, 
the swashplate mechanics (rotating and non-rotating 
components), pitch links, pitch horns, the rotor shaft 
and the hub (Fig. 1). The rotor blades are modeled 
as elastic beams undergoing coupled flap, lag and tor- 
sion deformation similar to the finite element meth- 
ods used in classical rotorcraft analyses. 

A third analytical model of the SASIP tiltrotor has 
been developed using a classical rotorcraft analysis 
known as UMARC/G, and is based on the UMARC 
(Ref. 5) comprehensive rotor code. This analysis 
does not have the capability to model complex joints 
and extreme nonlinear behavior as do the multibody 
codes, but is useful to serve as an analytical standard 
for some portions of the current study. 

In Ref. 6 the multibody models were correlated 
with experimental data from the SASIP model in 


hover, mainly concerning structural dynamics. The 
analysis focused on the investigation of the nonlinear 
behavior of the soft-inplane lead-lag hinge, on the 
interaction of the rotor motion with the control sys- 
tem, and on ground resonance stability issues. This 
paper extends that original work to the analysis of 
forward flight configurations, addressing whirl flutter 
and load prediction issues. 

The paper approaches the problem by directly 
comparing the results of the two multibody analysis 
codes, mainly focusing on those that are also available 
from the experimental campaign described in Ref. 7 
and on a selection of test cases that help speeding up 
model correlation, e.g. modal analysis of the rotor in 
vacuo. 

KEY RESULTS 

Several experiments have been conducted using the 
SASIP rotor system and its subcomponents (Ref. 7). 
The following list presents an overview of tests that 
are included in the analytical comparisons presented 
in this paper: 

1. Single blade cantilevered outside lag hinge: elas- 
tic mode comparison. 

2. Single blade mounted on hub in pendulum con- 
figuration, flap and lag hinge with lag spring and 
lag damper included: elastic mode comparison. 

3. Control system stiffness calibration: control stiff- 
ness and deadband comparison. 

4. Pitch-flap coupling and pitch-lag coupling cal- 
ibration: control system response comparison 

and analytical models. 

5. Hover run-up: fan plot comparisons. 

6. Hover performances: thrust, torque and blade 
angles as functions of collective at prescribed 
rpm; rpm sweep at prescribed collective setting. 

7. Wing structural analysis and hover stability sub- 
ject to ground resonance conditions: comparison 
of wing and rotor mode damping and frequen- 
cies. 

8. Airplane mode stability tests subject to 
whirlffutter conditions: comparison of wing and 
rotor mode damping and frequencies 

9. Conversion modes maneuver simulations: com- 
parison of dynamic load predictions. 

The experimental results for all nine of the listed 
comparisons have been obtained during several re- 
cent test campaigns; many of the data plots required 
for these comparisons were presented in Ref. 7. The 
analysis models have been refined through compari- 
son with the experimental data in sequence with the 
above list. 


Some of the key results of this investigation ob- 
tained to date are discussed in the following para- 
graphs. The multibody model developed for the 
three-bladed stiff-inplane wind-tunnel model is illus- 
trated in Fig. 1. The figure shows the rotor blades, 
pitch link, swashplate, and hydraulic control actua- 
tors which are attached to the pylon, and an elastic 
wing that is modeled using finite elements. For the 
SASIP four-bladed rotor system the rotor blades, hub 
joints, and most of the control system have been de- 
veloped. Figure 2 shows the DYMORE model of the 
SASIP setup. 

The wing and pylon model used for the three- 
bladed model is the same as for the SASIP rotor, so 
the wing/pylon submodel simply needed to be joined 
to the SASIP rotor model for the stability investi- 
gations. The aeromechanical stability of the coupled 
system will be investigated in future works; this work 
addresses the hover case. 

Cantilevered Blade Modal Analysis 

A complete review of the blade structural proper- 
ties has been conducted, to develop a reliable finite 
element model for subsequent analyses and accurate 
load recovery. Since both the multibody software 
codes allow for direct finite element modeling of ro- 
tor blades, this analysis permitted the determination 
of a reasonable trade-off between model accuracy and 
computational cost, leading to a reduced order model 
with respect to those used for detailed modal analysis 
in NASTRAN and UMARC/G. 

Table I reports the frequencies of a single blade, 
cantilevered right outside the lead-lag hinge. MBDyn 
results refer to a blade model made of 5 parabolic 3 
node C° beam elements (Ref. 8), with 11 structural 
nodes. DYMORE results refer to a 4 cubic beam 
FEM model, while in NASTRAN and UMARC/G 
25 beam elements were used. 

Hub-Mounted Blade Modal Analysis 

A comparison of elastic blade frequencies, for the 
condition of an isolated blade mounted to the hub 
(experimental setup shown in Fig. 3), is listed in Ta- 
ble II; the blade pitch is rigidly constrained. This 
analysis fully exploits the capabilities of multibody 
software in dealing with the exact kinematics of the 
articulated blade attachment and of the control sys- 
tem. Different configurations in terms of collective 
setting, imposed flap or lag angles, lead-lag hinge and 
pitch link stiffness have been addressed within a sin- 
gle model. The results indicate consistent capabili- 
ties of modeling the elastic blade and hinge dynamics 




Figure 3: Ground vibration test of an isolated blade 
with hub fixed. 



Figure 4: Mode shape comparisons for a non-rotating 
coupled flap-lag-torsion mode at approximately 64 
Hz. 




Figure 2: DYMORE model of the SASIP hub setup. 


Table I: Frequencies of the elastic modes for a single cantilevered blade. 


Mode 
No. Type 

Experiment 

Frequency, Hz 
NASTRAN UMARC/G 

MBDyn 

DYMORE 

1 

FI 

10.71 

10.61 

10.61 

10.61 

10.45 

2 

LI 

29.20 

29.11 

29.46 

29.03 

29.28 

3 

F2 

47.76 

51.16 

48.99 

50.45 

50.55 

4 

T1/F3 

107.29 

107.23 

107.08 

107.43 

107.28 

5 

F3/T1 

119.07 

121.49 

119.86 

116.78 

116.27 

6 

T2 

141.83 

177.19 

140.19 

162.69 

166.67 


Table II: Frequencies of the elastic modes for a single blade on a fixed hub. 


Mode 

No. Type 

Experiment 

Frequency, Hz 
NASTRAN UMARC/G 

MBDyn 

DYMORE 

1 

FI 

0.00 

0.00 

0.08 

0.67 

0.11 

2 

LI 

6.46 

6.54 

6.43 

6.32 

6.51 

3 

F2 

21.70 

19.48 

20.06 

19.37 

19.44 

4 

F3/L2 

61.15 

63.12 

64.20 

62.43 

64.30 

5 

T1 

107.94 

107.44 

103.50 

106.58 

107.07 

6 

F4 

119.25 

91.25 

96.21 

88.11 

92.30 


Control System Stiffness 



Figure 5: Equivalent control system stiffness; exper- 
imental data refers to measurements on each blade. 

among the analyses and generally good agreement 
with the experimental results. Mode shape compar- 
isons of the three analyses for the fourth mode are 
shown as an example because this mode has signif- 
icant participation from flap, lag, and torsion. As 
shown, the agreement in the flap and lag deflections 
is excellent, but there is a prediction discrepancy in 
the torsion participation. This difference in torsion 
dynamics is currently under investigation. 

Control System Calibration 

Direct measurements showed that the control sys- 
tem may be regarded as rigid up to the rotating 
swashplate. Measurable compliance has been de- 
tected only between the blade pitch and the swash- 
plate; moreover, an appreciable deadband (±0.4 deg) 
results, possibly related to bearing and pitch link ball 
joints wear. From the multibody analysis standpoint, 
this has been modeled by concentrating the compli- 
ance and the deadband into the rod that represents 
the pitch link. The resulting analytical equivalent 
control system stiffness (blade root torsional moment 
per unit blade pitch) is strongly dependent on the col- 
lective setting (Fig. 5). Moreover, the experimental 
data apparently show that there is a strong depen- 
dency on the direction the control system is loaded; 
this can be partly related to a significant presence of 
deadband. 

Control System Couplings 

The nonlinear modeling capability of the multi- 
body codes is highlighted in Fig. 6, which compares 
the experimental and predicted pitch-flap coupling 
response of the rotor system as a function of collec- 
tive, at zero flap and lag angles. The comparison 



Figure 6: Pitch-flap (S 3 ) and pitch-lag (64) coupling 
comparison between experiment and analysis. 

shows good trends, but a slight difference in magni- 
tude. The control system model is currently being re- 
fined to produce a better comparison with the exper- 
imental data. It is worth noting that the kinematic 
couplings also depend on the reference configuration 
about which they are computed; the only alternative 
to multibody exact kinematics, with less than ideal 
accuracy, is represented by pitch, flap and lag tabu- 
lation, and 3 parameter table lookup. 

Hover Run-Up 

Figure 7 shows the dependency on the collective 
of the frequencies of the blade mounted on the hub, 
without aerodynamics, when the deformability of the 
pitch link is accounted for. Both the non-rotating and 
the 100% hover angular velocity cases are considered. 
While the flap and the lead-lag rigid modes appear 
to be very little affected by the collective setting, the 
higher modes show more pronounced changes. Modes 
labeled as F3/L2 and F4 change their nature when 
subjected to the centrifugal effect; their shape is now 
dominated by the torsional deformation, and their 
denomination rather reflects the prevailing participa- 
tion of beam and chord bending. 

A plot of the regressive rotor lag mode frequencies 
as function of rotor speed is shown in Fig. 8. The plot 
shows a difference between the predicted and experi- 
mental results, although agreement for the nonrotat- 
ing condition is good. The lag hinge of the experi- 
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Figure 7: Dependency of isolated rotor in vacuo fre- 
quencies on collective settings; experimental results 
(thick hollow circles) refer to rigid pitch link; nu- 
merical results refer to deformable pitch link; thick 
squares represent DYMORE results at 10 deg collec- 
tive. 


experiment O 

MBDyn - tuned stiffness w/modal participation • 

DYMORE - crossover stiffness 

MBDyn - crossover stiffness ■ 

Figure 8: Regressive lead-lag mode frequency varia- 
tion with rotor speed. 

mental system is a complicated mechanism, and the 
results obtained thus far represent a simple, constant 
stiffness equivalent spring hinge joint model. Accord- 
ing to simple rigid blade theory, the regressive and the 
progressive lead-lag frequences based on root stiffness 
and rotation speed is 

V 

No unique lead-lag root stiffness could be found 
that matches the experimental results over the en- 
tire range of rotor rpm; Figure 8 shows the experi- 
mental results compared to the best fit obtained by 
calibrating the lead-lag root spring at the cross-over 
frequency, compared to those obtainable by calibrat- 
ing the spring at each rpm. A detailed modal anal- 
ysis showed a slight participation of blade chordwise 
bending in the mostly rigid lead-lag mode, resulting 
in a frequency reduction at 0 RPM of about 3.75%; 
this required an increase in the spring stiffness by 
about 7.5%. It is not yet clear what mechanism, if 
any, relates the lead-lag spring to the rotation speed. 

A complete fan plot of the rotor at 10° collective is 
presented in Fig. 9; some sensitivity of the rotating 
modes of the rotor to the collective appears at the 
low angles that are typical of hover; a detailed rotat- 
ing modal analysis at the higher collective settings of 
forward flight in a range of angular velocities will be 
the object of future analysis. The fundamental tor- 
sional mode is strongly dependent on the stiffness of 
the control system; the nominal frequency of about 
106 Hz with rigid control system shown in Table II 
drops to less than 100 Hz when the pitch link is mod- 
eled as a deformable rod, according to the calibra- 
tion described in Figure 5. The strong dependency of 
the control system stiffness on the collective setting 
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Figure 9: Fan plot of the isolated rotor in vacuo at 
10° collective; lines: DYMORE, symbols: MBDyn. 

is highlighted by the change in control/torsion fre- 
quency with the blade collective shown in Figure 7; 
note how the first torsion/control mode and the mode 
labeled as F4 cross at low and high collective. 

Hover Performances 

The hover performances of the SASIP rotor are de- 
tailed in Figs. 10-19. Experimental results are also 
shown, whenever available. The figures clearly show 
how the blade deformability highly affects the rotor 
performances. This is mostly related to blade torsion 
as a consequence of the offsets of the normal stress 
center, of the center of gravity and of the shear cen- 
ter of the blade sections. In fact, it has been verified 
that by artificially increasing the stiffness of the blade 
model, the curves shift towards those labeled as “rigid 
blade” . In the multibody analysis, simple uniform in- 
flow was used without any empirical correction, and 
a standard 2% total tip loss was considered. The 
thrust was not directly measured, since the WRATS 
rotor does not have an internal balance; the values in 
the figures refer to measures of the wing root beam 



Figure 10: Thrust vs. input collective. 

bending moment by means of strain gages, empiri- 
cally corrected to account for wing download. It is 
estimated that wing download amounts to 10% of the 
nominal thrust. 

Fig. 15 presents the actual blade pitch angle com- 
pared to the input collective. The differences between 
input and output angles mainly result from the con- 
tributions of the flap and lag angles through the con- 
trol system couplings. While for low input collective 
the flap and the lag effects are comparable in magni- 
tude but opposed in sign, and thus nearly cancel, for 
high input collective they add, resulting in an appre- 
ciable reduction in actual blade pitch. 

Figures 16 and 18 have been corrected by arbitrar- 
ily shifting the cone angle, since there were uncer- 
tainties on the reference value; Figures 17 and 19 
have been corrected as well by arbitrarily shifting the 
lead-lag angle, since the angle at rest was used as ref- 
erence, but there is no guarantee that it corresponds 
to the nominal blade rest azimuthal orientation. It 
is reasonable to assume that the rest position should 
be as close as possible to the more critical operating 
condition from a structural design standpoint, since 
the main reason for developing soft-inplane tiltrotors 
is to reduce the hub loads that are typical of stiff- 
inplane tiltrotors. This can be achieved by carefully 
tuning the lead-lag hinge chordwise position and the 
lead-lag spring orientation. 

Figs. 18 and 19, reported as a cross-check, show 
the expected rectilinear behavior that is typical of 
articulated rotors. 

Wing Structural Model and Hover Stability 

The wing model, in the case of hover and ground 
resonance analysis, in absence of significant sources 
of nonlinear behavior, does not present any peculiar 
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Torque vs. Collective; 875 rpm 



Figure 11: Torque vs. input collective. 


Thrust vs. Torque; 875 rpm 



Figure 12: Thrust vs. torque. 


Induced Velocity vs. Collective; 875 rpm 



Figure 13: Induced velocity vs. input collective. 


Thrust vs. Induced Velocity; 875 rpm 



Figure 14: Thrust vs. induced velocity. 


Blade Pitch vs. Collective; 875 rpm 



Figure 15: Blade pitch vs. input collective. 


Blade Cone vs. Collective; 875 rpm 



Figure 16: Blade cone vs. input collective. 
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Blade Lag vs. Collective; 875 rpm 



Figure 17: Blade lead-lag vs. input collective. 




Figure 20: Lead-lag modes vs. collective in hover 
(875 rpm). 
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Figure 18: Thrust vs. cone. 


Torque vs. Lag; 875 rpm 



Figure 19: Torque vs. lead-lag angle (positive when 
lagging). 


difficulties, as soon as the interaction of the rotor 
wake with the wing is not a concern. It is essential 
that the wing behavior, in terms of frequency and 
modal shapes at the interface with the rotor, is accu- 
rately modeled. The multibody formalism allows two 
different approaches to this problem: (a) the direct 
finite element modeling of the wing structure, and 
(b) a modal synthesis of the wing behavior. Within 
MBDyn both approaches have been considered, while 
with DYMORE the attention has been focused on the 
direct FEM modeling. 

The hover stability of the SASIP model is a sig- 
nificant issue. In fact, the SASIP experimental sys- 
tem, by design, is prone to ground resonance, because 
the soft-inplane rotor is mounted on a deformable 
wing whose fundamental frequencies are right below 
the non-rotating regressive lead-lag frequency, and 
the damping of the support, significantly in hover, 
is inherently delegated to the structure. At the same 
time, the lead-lag motion of the blades is damped at 
a nominal damping ration of 10-12%. 

Figure 21 shows the lead-lag angle rate that re- 
sults from a transient analysis of the entire model, 
with and without aerodynamics, during a linear RPM 
sweep at a 100 RPM/s rate. It is compared to the 
plot of some significant system frequencies in the 
fixed frame. The time and the RPM abscissae of 
the two subplots are proportional; this highlights 
interesting behavior patterns. At about 100 -j- 200 
RPM there seems to be some interaction between 
the regressive lead-lag, the progressive flap and the 
wing beam/chord frequencies. At 510 RPM the rotor 
speed crosses the lead-lag frequency, resulting in sig- 
nificant resonant response when the aerodynamics are 
considered. Finally, above 800 RPM the regressive 
lead-lag frequency crosses again the wing beam/chord 


Thrust vs. Cone; 875 rpm 



10 


PU 





Regressive lag ■■■•■ 




Wins beam/chord ■■ ■■ 







; II 


1 ... 



....a 




■ 

\t- 

"■•*1 

• 

■ 

•••••'-•« i • 



0 200 400 600 800 1000 

Rotor Speed [RPM] 



Figure 21: Lead-lag angular velocity during RPM sweep at 10 deg collective. 


mode; this results in ground resonance when aero- 
dynamics are neglected, because no wing structural 
damping was considered. The aerodynamics seem to 
provide enough damping to the system to prevent dy- 
namic instabilities. 

Airplane Mode Stability 

The new four-bladed, soft-inplane rotor system, 
oriented in airplane mode for high-speed wind- 
tunnel testing, is shown in Fig. 22 mounted on the 
WRATS model in the NASA Langley Transonic Dy- 
namics Tunnel (TDT). The basic dynamics of the 
wing/pylon/rotor system shifts substantially with 
conversion to airplane mode, as the mass offset of 
the pylon/rotor moves from above to forward of the 
elastic axis, and thus creates a significant coupling 
between the wing beam and torsion modes and the 
rotor lag mode. The wing chord mode becomes pre- 
dominantly isolated from these modes in the airplane 
configuration. 

For airplane-mode aeroelastic stability testing, the 
rotor system is normally operated windmilling (un- 
powered and disconnected from the drive system), 
with the collective blade pitch used to adjust the ro- 
tor speed, and there is near-zero torque at the rotor 
shaft. This represents the most conservative man- 
ner to test the stability of the system (no damping 
from the drive system). Under windmilling operation, 
damping of the key mode associated with system sta- 



Figure 22: WRATS SASIP in forward flight configu- 
ration in the Langley Research Center Transonic Dy- 
namics Tunnel (TDT). 




Figure 23: Wing beam bending mode damping com- 
parison between original stiff-inplane WRATS model 
and SASIP soft-inplane model. 

bility (the wing beam mode) was determined to be 
significantly less for the new four-bladed soft-inplane 
hub than for the three-bladed stiff-inplane (baseline) 
system, as shown in Figure 23 (from Ref. 7). Damp- 
ing of the wing beam mode was generally less than 
1.0% in windmilling flight for all the soft-inplane con- 
figurations considered (on-downstop (D/S), off-D/S; 
0.57/rev dampers, 0.63/rev dampers; 550, 742, and 
888 RPM rotor speeds). 

In powered-mode (200 in-lb torque maintained, 
~22.6 Nm) the system damping and the stability 
boundary both increased significantly as illustrated 
in Figure 24 (again from Ref. 7; note that the on- 
D/S configuration is shown rather than the off-D/S 
as used in Figure 23 because of low damping associ- 
ated with the off-D/S case). 

Figures 25-28 show the frequency and the damping 
of the wing beam mode at different rotor speeds in 
off- and on-D/S configuration. The numerical values 
and trends are compared to the results from the wind 
tunnel campaign described in Ref. 7. 

Conversion Maneuver Simulation 

The wing-pylon attachment setup of the WRATS 
SASIP model does not allow the simulation of the 
conversion maneuver; different conversion angles can 
be obtained in the wind tunnel by changing the length 
and the stiffness of the rod that models the conver- 
sion actuator. The stiffness of the connection ele- 



Figure 24: Wing beam bending mode damping com- 
parison between powered and windmill configurations 
of the WRATS SASIP soft-inplane model (solid sym- 
bols are MBDyn analysis). 

ment is usually modified also to reproduce the conver- 
sion mechanism stiffness in the on- and off-downstop 
configurations. However, the availability of validated 
multibody models of the WRATS SASIP system al- 
lows to simulate rather complex manuevers, includ- 
ing those, like the conversion, that imply finite, and 
definitely large, rotations and changes in system con- 
figuration. The problem has been addressed by the 
multibody analysis, considering a complete conver- 
sion performed at a constant rate of 3 deg/s at differ- 
ent airspeeds, to obtain an estimate of the structural 
loads within the conversion corridor. 

COMPUTATIONAL ASPECTS 

The two multibody analyses considered in this 
work present common aspects and some differences. 
The approach to the analysis of rotorcraft systems 
is basically analogous; it is partially different from 
usual comprehensive rotorcraft analysis, since it is es- 
sentially based on performing a virtual experiment : a 
detailed, nonlinear model of the rotorcraft system is 
analyzed by repeatedly performing time integration 
of Initial Value Problems (IVP), and synthetic infor- 
mation, like thrust or torque levels, blade angles at 
trim are obtained by averaging the results of the sys- 
tem response, or by letting it reach a steady solution, 
if any. Modal analysis follows different approaches. 




WRATS SIP T565; OFF-DS 550 RPM 
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WRATS SIP T565, OFF-DS 550 RPM 



Figure 25: Wing beam bending mode frequency and 
damping comparison between experimental and nu- 
merical solutions at 550 RPM, off-D/S in windmill 
conditions 


WRATS SIP T565; ON-DS 560 RPM 




Velocity, knts 

Figure 26: Wing beam bending mode frequency and 
damping comparison between experimental and nu- 
merical solutions at 550 RPM, on-D/S in windmill 
conditions 
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Figure 27: Wing beam bending mode frequency and 
damping comparison between experimental and nu- 
merical solutions at 742 RPM, off-D/S in windmill 
conditions (solid symbols are MBDyn analysis). 


Figure 28: Wing beam bending mode frequency and 
damping comparison between experimental and nu- 
merical solutions at 742 RPM, on-D/S in windmill 
conditions (solid symbols are MBDyn analysis). 


DYMORE uses an implicit matrix method that ex- 
ploits the properties of the Arnoldi’s algorithm, a 
subspace method, which needs only a matrix- vector 
multiplication to extract the highest modulus eigen- 
values of a discrete system, e.g. those that lie outside 
the unit circle in the complex plane (Ref. 9). 

MBDyn, instead, uses the Proper Orthogonal De- 
composition (POD) to extract a set of basis functions, 
called Proper Orthogonal Modes (POM), from the re- 
sults of the numerical simulations. They are sub- 
sequently used in a Galerkin projection that yields 
low-dimensional dynamical models. The POMs are 
a minimal set of output signals that can be used to 
identify the dominant eigenvalues of the transition 
matrix (Refs. 10, 11). The numerical integration of 
the underlying IVP is performed by means of an orig- 
inal implicit multistep integration scheme, that guar- 
antees second-order accuracy with tunable algorith- 
mic dissipation. 

Either of the two techniques is required because 
the detailed tiltrotor multibody model can be quite 
large; for instance, the complete MBDyn model en- 
tails about 800 unknowns, and models with as much 
as 2000 unknowns can be quite common. The vir- 
tual experiment approach requires the execution of 
long runs, so its feasibility heavily relies on the effi- 
ciency of the software. The codes used in this work 
proved to be able to perform the required computa- 
tions, with rather realistic and detailed models, in 
reasonable times, as shown in Table III, along with 
other interesting figures. 

The analysis of specific configurations and flight 
regimes may require the determination of non-trivial 
trim points. This work required the analysis of two 
significant cases: the windmill and the powered con- 
dition. 

Windmill: this operational regime is defined by 

the collective setting that allows to maintain constant 
RPM with no power; this collective is unknown, so 
to obtain the trim point, a low-gain controller that 
changes the collective based on the integral of the 
RPM error has been implemented. A simplified form 
of the mast equilibrium equation is 

jfl = -C. 

Its linearization yields 

r)C r)G 

jn + ^(n-n 0 ) + ^(e-e 0 ) = -c : (i) 

where the sensitivities of the torque C to the angular 
velocity SI and the collective 9 are both positive. The 
controller equation is 


which, combined with Eq. (1), normalized by the in- 
ertia J, yields 



-C /n -C/g 

w 

G 0 


{ n } + { C/oDo + C/A } 


It is stable provided positive gains G are considered; 
the poles are 

C/n T Jc^-AGC,g 
S ~ 2 

The dominating one is that with the minus sign be- 
fore the square root; for G small enough, it can be 
approximated by 


C/g 

s ^-^Ag-, 

C /u 


this eigenvalue dictates the time constant. Note that 
C/g depends on the square of the airstream speed, 
while C/ci is virtually independent from it; as a conse- 
quence, the time constant r = —1/s for a given gain 
G decreases as the airstream increases. Time con- 
stants above 1 s have been used, to clearly separate 
the dynamics of the governor from the wing dynam- 
ics, in general above 5 Hz. The controller drives the 
collective to the value that corresponds to zero thrust 
in a relatively short time (10-^30 revolutions, depend- 
ing on the airstream speed through the sensitivity of 
the aerodynamic torque to collective changes), while 
preserving the windmill degree of freedom of the hub. 
Powered: this operational regime has been simu- 

lated in two different manners. First, a setup anal- 
ogous to that of the windmill regime has been used, 
with the desired torque applied between the shaft and 
the pylon. As a consequence, the controller drives the 
collective to a value that provides an aerodynamic 
torque equal to the applied one, while preserving the 
RPM. This case represents a powered trim point with 
the windmill degree of freedom still in place. Another 
case has been analyzed, consisting in imposing the 
RPM, thus eliminating the windmill degree of free- 
dom, with the collective setting found above. This 
latter case corresponds to an infinite gain governor. 
As a general trend, it is observed that the presence 
of some torque has a slightly stabilizing effect, which 
increases further by imposing the RPM. Apparently, 
of the two powered cases ranging from powered trim 
with windmill degree of freedom to infinite gain gov- 
ernor, the latter is more representative of the WR ATS 
SASIP wind-tunnel model. 


e = G{n-n 0 ) 



Table III: Numerical Aspects 


Data 

DYMORE MBDyn 

Model Size 

~800 unknowns 

Beam Elements per Blade 

5 

Degrees of Freedom per Blade 

132 

Wing Modes 

5 

Wing Beam Elements 

— 

Time Step 

0.0014-0.0005 s 

Time Steps per Revolution 

804-160 

Real/Simulated Time Ratio 

~45:1 

Computer 

Athlon 2.2GHz 


CONCLUDING REMARKS 

This study shows that multibody codes can be suc- 
cessfully used to model complex mechanisms in rotor- 
craft analysis and that this capability can improve 
predictions of dynamic behavior. How this improved 
modeling capability can influence load and stability 
predictions is an anticipated conclusion which has not 
yet been resolved. Key points that highlight the im- 
portance of exploiting the capabilities of multibody 
analysis software are: 

• hub and blade root kinematics and dynamics 

• control system kinematics and structural dynam- 
ics 

• finite element and composite-ready blade struc- 
tural modeling 

• versatility in structural modeling of unusual 
components 

• essential rotor aerodynamics capabilities 

The two multibody analysis software show analo- 
gous capabilities when applied to rotorcraft modeling. 
They are both able to capture the essential aspects of 
the whirl flutter stability of the WRATS SASIP wind- 
tunnel model, although the relatively inaccurate wing 
aerodynamic models available, and some uncertainty 
on the level of structural damping, did not allow, in 
some specific case, to obtain consistent quantitative 
whirl flutter predictions. The availability of two inde- 
pendently developed and validated multibody mod- 
els of the WRATS SASIP represents an opportunity 
for future investigation of the aeromechanics of soft- 
inplane tiltrotor systems. 
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